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0^ ■ Abstract 

"^ I A critical point of second order, belonging to the universality class 

o ' of the 3d Ising model, has recently been advocated as a strong can- 

ji. . didate for the critical behaviour (at high temperatures) of QCD with 

D I non-zero quark masses. The implications of this conjecture are inves- 

i-^ ' tigated in the multiparticle environment of high-energy collisions. A 

universal intermittency pattern of pion-density fluctuations is found, 
k> i at the critical point, and its association to the critical exponents is 

5h ' discussed. A Monte Carlo simulation of critical events, in heavy-ion 

collisions, reveals the detailed structure of these fluctuations, suggest- 
ing a framework of (event-by-event) measurements in which the critical 
theory of QCD may become falsifiable. 
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The nature of the QCD critical behaviour, associated with chiral phase 
transition, has gradually been better understood, during the last few years, 
in terms of the number of flavours and the primordial mass spectrum of 
the light quarks. In an approximate world of two flavours {u, d) with zero 
quark masses (m„ k, nid k, 0), the high-temperature QCD phase transition 
(Tc ^ ISOMeV) is of second order [1-3]; it belongs to the 0(4) universality 
class of a ?>d Heisenberg magnet [^ and the order parameter is associated 
with pion and sigma field condensates = (vr, a). In a more realistic world 
with three flavours (m, d, s) and a finite mass {ms) of the strange quark, the 
system develops a tricritical point (m* , T^) at which the transition changes 
from first {iris < 1^*3) to second order {ms > m*) @, ^. In the real world, 
however, with non-zero quark masses {mu,md 7^ 0) the second-order branch 
disappears from the phase diagram: the corresponding transition becomes a 
smooth crossover [Q and the tricritical point is transformed into an endpoint 
of the first-order phase transition line (critical point of second order). It 
has also been argued [^, |^ that the behaviour of the system near the QCD 
critical point {Tc, fj,c) belongs to the universahty class of the Ising model in 
3d (gas-liquid transition); in this case, the order parameter is associated 
with the cr- field condensate alone, a ~ {'4''4')- The above over-all picture 
of critical QCD has been advocated recently by a number of authors [3- 
8] and despite the fact that one cannot (yet) verify these ideas completely 
in a detailed ab initio study on the lattice, owing to the involvement of 
a non-zero chemical potential /ic at the critical point, this picture remains 
attractive, as it is based on rather general and universal properties of strong 
interactions, near criticality. In this Letter we claim that it can also be 
falsifiable, because it may lead to quantitative and sharp predictions for 
pion-density fluctuations, near the critical point. Such critical fluctuations 
can, in principle, be measured in relativistic heavy-ion collisions, by selecting 
rare events with strong and universal intermittency pattern in momentum 
space. 

To this end, consider the 3d effective action rc[o"], which effectively de- 
scribes the system of QCD at the critical point (T = Tc, fi = ^c) '■ 



Tc[a] = T^'Jd'x [^{Vaf + GT^{T-'af^' 



(1) 



In eq.(|l]) the macroscopic field a has a dimension a ~ (length) ^, 5 is the 
isothermal critical exponent, and G, a dimensionless coupling, specifies the 

1 



critical equation of state: -^ ~ G(f)^. Both parameters (G, S) are universal; in 
the 3d Ising class one may fix the exponent at the mean field value 6^5, be- 
cause of the smallness of the anomalous dimension []TU|. The coupling G, for 



the same universality class, has been estimated recently, in non-perturbative 
studies of the critical equation of state, and its actual value is fixed, to a good 
approximation, in the range G ~ 1.5-2.5 [|10|. Hence the only free parameter 
in the effective theory (|l|) is essentially the critical temperature T^, which 
determines the length scale (in heavy-ion physics, /3c = T~^ ~ 1-2 fm). 

In order to test the theory (|l]) in the environment of a high-energy colli- 
sion, we have to first adapt the longitudinal geometry to the rapidity space 
(^) by introducing the proper time-scale (r) of the collision [A + A): 



L//1 J 




2 iy^aY + Gcipy 



(2) 



where Ca = ^- The effective field (t(^,x_l), being macroscopic (classi- 
cal), is naturally associated with the density of cr-particles in the 3(i space 
(^,x_l): (J^ = A^o- (volume)"^ ~ (length)"^. Technically, this interpretation 
follows from the requirement that the density matrix, associated with the 
partition function Z = J[Da]e~^'' of the critical theory (|^), be diagonal in 
the coherent-state representation of the particles associated with the a-field 
[ p^ . As a consequence, one finds for the average multiplicity (iVo-) the form: 
{N„) = Z~^ j\Da]e~^''{j d'^x^d^ cr^(^,xj_)), suggesting that the appropriate 
order parameter for the study of the density fluctuations of the isoscalar 
particles at the critical point, is o"^(^,afx) [ll2[ . 



In the cylindrical geometry of the collision, the theory (|^) can be solved 
by projecting out the effective action onto rapidity and transverse space. 
Following the results in ref. ]13|, the critical fluctuations of the order pa- 
rameter in these two spaces {Id and 2d) are developed within clusters of size 
5c = "" ^^^ and r_Lc = '^ ^^^ respectively {R± is the transverse radius of the 
entire system at the critical point). These fluctuations are properly described 
(with the expected fractal dimension) by the contribution to the partition 
function of instanton-like configurations, which, within these clusters, are 
slowly varying functions and contribute equally to the effective potential and 
to the derivative terms in eq. (0) [l^. In this framework, the canonical 
partition function, at T = Tc, for scalar particles (cr) in a cylindrical volume 



V < 7!'6cr'j_^, is written: 



Z{N,,V,T,) = N;'/'exp 



¥-' 



(3) 



where K = (31./2GC2. 

The average multiphcity {N^) of scalars, in each Sd cluster, as a function 
of the volume V , now follows from @: 

Equation (^ shows that the critical clusters are local, cylindrical fractals 
with a minimal volume-scale V = Vo where self-similarity breaks down. The 
fractal dimensions in rapidity and transverse space are dp^i = | and dF,2 = |, 
reflecting the actual value of the isothermal critical exponent 5 ~ 5 and 
also the choice of the order parameter a^ {dp^i = |^, dF,2 = \^i ) ■ The 
corresponding density- density correlation function obeys for each cluster a 
characteristic power-law: 

w,..-j,(o.o)>^(-||-)'%-ffi|)|r"»i|i-" (5) 

Since the size of the critical clusters specifies also the correlation length of 

— 1/9 

the finite system |13|, both in rapidity {6c = 2t^ (§) ) ^^^ transverse 



space, (r_L c = ^^{^)~^^'^R±) we may safely assume that in the global volume 
Vg = 7rR'j_A (A is the total rapidity available in the experiment) the critical 
clusters interact weakly, so that they form an ideal classical gas with a Poisson 
multiplicity distribution [|14|. The average multiplicity of clusters, in this 



,2 



approximation, is {Nd) ~ (7^) [t] ^^d, using eq. (§), an upper bound 
of the excess multiplicity of pions produced when the system passes through 

the critical point, can be extracted (A^^ ^ 2N^) ^ < ^ (^7) (f )'^ ^^ 
On the other hand, the local fractals, associated with the critical clusters 
finally build-up a global fractal with the same dimension [|l^. In fact, in 
the transverse momentum space, the fractal F2 is even stronger, since the 
corresponding Fourier transform of the power-law (H) gives dF,2 = | ('^f,2 = 
2 — dp^2)- On the basis of these structures, the theory predicts a universal 



pattern of strong, multidimensional intermittency (in rapidity, transverse 
momentum and 3d momentum space) of the form: F^^^ ~ M~3~, F(^) ~ 



2(p-i) ^^o^ . , ,Qx 5(p-i) 



[M'^)^^, Fj(3) ~ {M^)^r- {p = 2,3,.) where Fp is a scaled factorial 
moment of order p as a function of M , the number of one-dimensional cells 
in momentum space. Other power-laws coming from the Id projections of 
the fractal F2 or from the Cartesian product of two fractals in the cylindrical 
geometry of the system, can also be written [1^. The predicted intermittency 



power-laws describe the density fluctuations of scalars produced at the critical 
point quantitatively, relating these fluctuations to the critical exponents of 
the 3d universality class. 

We have performed a Monte Carlo simulation of the critical events in 
Pb-(-Pb collisions in the SPS (A ^ 6), corresponding to the choice of the 
parameters: T^ = 120 MeV, R^ = 30 fm, Tc = 25 fm. We have also fixed 



the coupling G at the value G = 2 ||10[. In momentum space, the isoscalar 



particles (cr) are grouped in critical clusters (cyhndrical) with appropriate 
power-law correlations, dictated by the Fourier transform of eq. @ in trans- 
verse space: (p(^,p_l)p(0, 0)) ~ \^\~^^^\PcP±\~'^^^ ■ The distribution of the 
clusters in rapidity (projected) corresponds to an ideal Id gas, whereas in 
transverse momentum space the following extra constraints are imposed: (a) 
the clusters (their centers) are exponentially distributed according to the 
temperature Tc, dric ~ exp(— ^y^) and (b) the total transverse momentum 

carried by the clusters is zero, J2iP±c — 0- Moreover the size of the criti- 
cal clusters in transverse momentum space, as determined by the minimum 
length scale in transverse configuration space, is given by the expression 

|Ap_l| ~-2^(-;^) (2) (^) ^^"^ f°^ t^^ above set of parameters we 
estimate |Apj_| ^ 70 MeV. With these ingredients and the overall constraints 
on the formation of sigmas at T = Tc, a typical, simulated event is obtained 
with characteristics shown in Figs. 1-3: the overall picture is consistent with 
a 3d cylindrical fractal in momentum space (Fig. 1), which is decomposed 
into a 2d self-similar structure in transverse momentum space (Fig. 2) and a 
Id self-similar structure in rapidity space (Fig. 3). The multidimensional in- 
termittency pattern clearly reveals a linear spectrum of indices in each case, 
consistent with the critical exponents of the theory |13 . 



In the process of a phase transition in thermal equilibrium, the fluctuating 
system is very likely to get frozen immediately after crossing the critical 



point of second order {^c,Tc) [0; therefore, the fluctuations predicted at 
T = Tc may enter the horizon of observation (in the corresponding events) 
without any distortion (Tj ^ T^) . It remains to guarantee that the decay of 
0" may proceed near the vrvr threshold, rather rapidly transferring the critical 
fluctuations found at T = T^ automatically to the vr" spectra (a -^ tt+vt") 
(Figs. 1-3). The isoscalar a has a very small mass at T = Tj; therefore, 
the system of sigmas remains frozen and stable until the mass m„, as a 
function of temperature, rises to the vrvr threshold value, ma = 2mT^. If the 
corresponding time scale Tth is sufficiently long (in particular events), the 
expansion rate |^| < 3:^, for r > Tth, becomes very small with respect to 
the decay rate of a (even near the tttt threshold), also because the isoscalar 
is very strongly coupled to the vrvr system. Under these circumstances, the 
sigmas with rucr ~ 2m.,r have enough time to decay, before their mass increases 
further, leading to charged pions (vr^) with the same fluctuation pattern in 
momentum space. We have verified numerically this decay mechanism, near 
the TTTT threshold, in the case of the typical event shown in Figs. 1-3. We 
have used for the parameters of the cr- meson in vacuum (T = 0) the average 
values m(,°) « T^^^ ^ 0.8 GeV of the /o(400 1200) vrvr resonance 0. It 



turns out that, if Tth ^ 10 fm, the decay a ^ nn proceeds near vrvr threshold 
and therefore the intermittency pattern found at T = Tc for sigmas (Figs. 
1-3) is also a prediction for the density fiuctuations of the abundant charged 
pions (tt") associated with the critical point. 

In the framework of heavy-ion physics, we have derived a pattern of lo- 
cal dynamical fiuctuations of the classical a-field associated with the critical 
point in QCD with non-zero quark masses. This critical point belongs to 
the 3d Ising universality class, but the same pattern is valid in any univer- 
sality class associated with a 3d scalar theory. The reason is that, in these 



theories, the anomalous dimension 77 is very small |]10|, leading to the same 
(mean field) value of the exponent: 5 ~ 5. Hence a similar pattern of fiuctu- 
ations is expected, for instance, in the 0(4) 3d theory [Q, which describes 
the chiral QCD phase transition in the limit of zero quark masses 0. At 
the level of observed pions, however, these theories are expected to lead to 
different fluctuation patterns since the in-medium evolution (and decay) of 
scalars varies from one universality class to the other. In the case of real 
QCD with non-zero quark masses, we have argued that in a particular class 
of events, in which the in-medium evolution of a is rather slow, Tth ^ ^t 



c; 



the fluctuation pattern found at T = Tc (Figs. 1-3) is inherited by the 
charged decay products of a near the tttt threshold (vr") and manifests itself 
as a clear signature of the quark-hadron phase transition in nuclear collisions. 

We thank A.L.S. Angelis and CG. Papadopoulos for their comments and 
suggestions on the present work. 
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Figure Captions 



Figure 1: (a) The distribution of the a-particles of a typical critical event 
in 3d momentum space. The multiplicity of a's is 138. (b) The corresponding 
{3d) factorial moments of order p = 2,3,4 in a log-log plot. A linear fit 
(dashed line) indicating the intermittency pattern for this event is included. 

Figure 2: (a) The distribution of the a-particles in transverse-momentum 
space for the event of Fig. 1. (b) The corresponding {2d) factorial moments 
of order p = 2, 3, 4 in a log-log plot. A linear fit (dashed line) indicating the 
resulting intermittency pattern is also shown. 

Figure 3: (a) The distribution of the a-particles in rapidity for the event 
of Fig. 1. (b) The rapidity factorial moments of order p = 2, 3, 4 in a log- 
log plot together with a linear fit (dashed line) indicating the corresponding 
intermittency pattern. 
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